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1. Introduction 

We introduce the class of V- variable fractals. This solves two major problems previously 
restricting further applications of Iterated Function System (IFS) generated fractals, 
the absence of fine control on local variabihty and the absence of a fast algorithm for 
computing and accurately sampling standard random fractals. 

The integer parameter V controls the number of distinct shapes or forms at each 
level of magnification (Figures 4, 6). The case V = 1 includes standard deterministic 
fractals (Figures 1, 2, also the fern and lettuce in Figure 4) generated by a single IFS 
and homogeneous random fractals (c.f. Hambly 1992, 2000 and Stenflo 2001). Large 
V allows rapid approximations to standard random fractals in a quantifiable manner, 
and the approximation is to not one, but to a potentially infinite sequence of correctly 
distributed examples. The construction of F-variable fractals is by means of a fast 
Markov Chain Monte Carlo type algorithm. In particular, one can now approximate 
stndard random fractals, together with their associated probability distribution, by 
means of a fast forward algorithm. 

A surprising but important fact is that each family of l^- variable fractals, together 
with its naturally associated probability distribution, forms a single superfractal 
generated by a single supcrlFS, operating not on points in the plane (for example) as 
for a standard IFS, but on ^-tuples of images. Furthermore, dimensions of V^-variable 
fractals are computable using products of random matrices and ideas from statistical 
mechanics. We implement a Monte Carlo method for this purpose. 

Since the mathematically natural notion of a ^-variable fractal solves two major 
problems previously restricting wider applications, we anticipate that y- variable fractals 
should lead to significant developments and applications of fractal models. For example, 
IFSs provide models for certain plants, leaves, and ferns, by virtue of the self-similarity 
which often occurs in branching structures in nature. But nature also exhibits 
randomness and variation from one level to the next — no two ferns are exactly alike, and 
the branching fronds become leaves at a smaller scale. V^-variable fractals allow for such 
randomness and variability across scales, while at the same time admitting a continuous 
dependence on parameters which facihtates geometrical modeUing. These factors allow 
us to make the hybrid biological models in Figures 4, 9. Because of underlying special 
code trees (Barnsley, Hutchinson and Stenflo 2003a) which provide the fundamental 
information-theoretic basis of V^-variable fractals, we speculate that when a l^-variable 
geometrical fractal model is found that has a good match to the geometry of a given 
plant, then there is a specific relationship between these code trees and the information 
stored in the genes of the plant. 

In this paper we describe the algorithm for generating V-variable fractals, explain 
the notion of a superfractal, and show how to compute the dimension of F-variable 
fractals in case a uniform open set condition applies. In order to make the ideas clearer, 
and hopefully facilitate the application of these notions to non-mathematical areas, we 
have avoided technicalities and illustrated the ideas by means of a number of model 
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examples, but the general constuction and results should be clear. 

In Barnsley, Hutchinson and Stenflo (2003a) we survey the classical properties 
of IFSs, develop the underlying theory of V-variable code trees and establish general 
existence and other properties for V^-variable fractals and superfractals. In Barnsley, 
Hutchinson and Stenflo (2003b) we prove the dimension results for l^-variable fractals, 
for which we here give an informal justification. 

2. Iterated Function Systems 

By way of background we first recall the concept of an IPS via the canonical example 
of the Sierpinski triangle S (approximated in the bottom right panel of Figure 1), 
which has been studied extensively (see Falconer 1990, 1997 and Hambly 1992, 2000) 
both mathematically and as a model for diffusion processes through disordered and 
highly porous material. The set 5* has three components, each of which is a scaled 
image of itself; each of these components has three sub-components, giving nine scaled 
images of at the next scale, and so on ad infinitum. A simple observation is that if 
fi, f2, fs are contractions of space by the factor | with fixed points given by the three 
vertices ^1,^2, A3 respectively of S, then the three major sub-components of S are 
fi{S)j2{S), fsiS) respectively and 

^ = /l(^)U/2(^)U/3(^). (1) 




Figure 1. Convergent or Backward Process. Beginning from any set (fish) Tq, 
iterates Ti = F{To), T2 = F{Ti),. . . converge to the Sierpinski Triangle S. Shown 
are iterates Tq, Ti, T2, T3, T4, Tg. 

The collection of maps F = (/i,/2,/3) is called an Iterated Function System (or 
IFS). For any set T one similarly defines F{T) = fi{T)Uf2{T)Ufs{T). It is not difficult 
to show that S is the unique compact set satisfying . Furthermore, beginning from 
any compact set Tq, and for k > 1 recursively defining = F{Tk-i), it follows that 
Tfc converges to S in the Hausdorff metric as ^ 00, independently of the initial set 
To (Figure 1). For this reason, S is called the fractal set attractor of the IFS F and 
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this approximation method is called the convergent or backward process, c.f. Hutchinson 
(1981). 

An alternative approach to generating S" is by a chaotic or forward process (Barnsley 
and Demko 1985), sometimes called the chaos game (Figure 2). Begin from any point Xq 
in the plane and recursively define Xk = fk{xk-i), where each fk is chosen independently 
and with equal probability from (/i,/2,/3)- With probability one the sequence of 
points {xk)k>o approaches and moves ergodically around, and increasingly closer to, 
the attractor S. For this reason F is called an iterated function system. If instead the 
fk are selected from {fi, f2, fs) with probabilities (^1,^2,^3) respectively, where each 
Pi > and pi + P2 + Ps = 1, then the same set S is determined by the sequence {xk)k>o, 
but now the points accumulate unevenly, and the resulting measure attractor can be 
thought of as a greyscale image on S, or probability distribution on S, or more precisely 
as a measure. In this case (/i, /2, f3',Pi,P2,P3) is called an IFS with weights. 




Figure 2. Chaotic or Forward Process. Beginning from any initial point and 
randomly and independently applying /i, /2 or f^ produces the Sierpinski Triangle 
as attractor with probability one. Shown are the first 10,000 and 100,000 points 
respectively. 

These ideas and results naturally extend to general families of contraction maps and 
probabilities. Even with a few affine or projective transformations, one can construct 
natural looking images (see the initial fern and lettuce in Figure 4). 

IFSs have been extended to study the notion of random fractals. See Falconer 
(1986), Graf (1987) and Mauldin and Williams (1986); also Hutchinson and Riischendorf 
(1998, 2000) where the idea of an IFS operating directly in the underlying probability 
spaces is used. 

3. Construction of V-variable fractals 

We now proceed to the construction of V^-variable fractals. This can be understood in 
the model situation of two IFSs F and G, with F as before and with G having the same 
fixed points Ai,A2,A3 but with contraction ratios | instead of |. We emphasise that 
the following construction is in no sense ad hoc, but is the natural chaotic or forward 
process for a superfractal whose members are V^-variable fractals as we see later. 



V -variable fractals and superfractals 



5 








GA 




F 


A ^ 


AA 


A 




AA 












F A ' 

/A 


G * ^ 

A 

\ A 


A \ 




F 

▲ ▲ 



Figure 3. Forward algorithm for generating families of ^-variable fractals. Shown 
are levels 1 (top), 2 and 3 in the construction of a potentially infinite sequence of 
5-tuples of 5-variable Sierpinski Triangles. For each buffer from level 2 onwards, 
F or G indicates which IFS was used, and the input arrows indicate the buffers to 
which this IFS was applied. 

One begins with arbitrary sets, one in each of V input buffers at level 1 (Figure 3 
where V = 5, colour coded to indicate the maps involved, and Fig. 6 where V = 2). A 
set in the first of V output buffers (level 2, Figure 3) is constructed as follows: choose 
an IFS F or G with the preassigned probabilities or respectively; then apply 
the chosen IFS to the content of three buffers chosen randomly and independently with 
uniform probability from the input buffers at level 1, allowing the possibility that the 
same buffer is chosen more than once (thus one is performing uniform sampling with 
replacement). The resulting set is placed in the first buffer at level 2. The content of 
each of the remaining buffers at level 2 is constructed similarly and independently of 
the others at level 2. These output buffers then become the input buffers for the next 
step and the process is repeated, obtaining the bottom row in Figure 3, and so on. 

The construction produces an arbitrarily long sequence of l^-tuples of approximate 
variable fractals associated to the pair of IFSs (F, G) and the probabihties (P^, P*^). 
The degree of approximation is soon, and thereafter remains, within screen resolution 
or machine tolerance (Figure 6). The empirically obtained distribution of V^- variable 
fractals over any infinite sequence of runs is the same with probability one and is the 
natural distribution as we explain later. The generalisation to the case of a family of 
IFSs (F^, . . . , F^) with associated probabilities (P^, . . . , P^) is straightforward; also 
sets can be replaced by greyscale images or more generally by coloured sets built up 
from primary colours if IFSs with weights are used (Figure 9). 
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Figure 4. 2-Variable Fractals. In the first "Squares and Shields" example, the 
number of components at each level of magnification is 1,4,16,64, . . . , but at each 
level there are at most two distinct shapes/forms up to affine transformations. The 
actual shapes depend on the level. In the second example of fractal breeding, fern 
and lettuce fractal parents are shown with four possible 2-variable offspring. The 
two IFSs used are those generating the fern and the lettuce. The associated 
superfractal is the family of all possible offspring together with the naturally 
associated probability distribution. 

The V -variability can be understood as follows. In Figure 3 the set in each buffer 
from level 2 onwards is composed of three component parts, each obtained from one 
of the V = 5 buffers at the previous level; at level 3 onwards sets are composed of 
9 smaller component parts each obtained from one of the V buffers two levels back; 
at level 4 onwards sets are composed of 27 smaller component parts each obtained 
from one of the V buffers three levels back; etc. Thus at each level of magnification 
there are at most V distinct component parts up to rescaling. In general, for a V- 
variable fractal, although the number of components grows exponentially with the level 
of magnification, the number of distinct shapes or forms is at most \^ up to a suitable 
class of transformations (e.g. rescalings, affine or projective maps) determined by the 
component maps in each of the IFSs {F^, . . . ,F^) (Figure 4 first panel. Figure 6). If 
parts of the V^-variable fractals overlap each other, the V^-variability is not so obvious, 
and fractal measures (greyscale or colour images) rather than fractal sets (black and 
white images) are then more natural to consider. 



4. Superfractals 

The limit probability distribution on the infinite family of V^-variable fractals obtained 
by the previous construction is independent of the experimental run with probability 
one as we explain later. As noted previously, an initially surprising but basic fact 
is that this family of \^-variable fractals and its probability distribution is a fractal 
in its own right, called a superfractal, and the construction process for the generated 
collection of V^-variable fractals turns out to be the forward or chaotic process for this 
superfractal, see also Barnsley, Hutchinson and Stenflo (2003a). For large V, l^- variable 



V -variable fractals and superfractals 



7 



fractals approximate standard random fractals and their naturally associated probability 
distributions in a quantifiable manner, providing another justification for the canonical 
nature of the construction, see Barnsley, Hutchinson and Stenfio (2003a, b). Although 
there was previously no useful forward algorithm for standard random fractals such as 
Brownian sheets, one can now rapidly generate correctly distributed families of such 
fractals to any specified degree of approximation by using the previously described fast 
forward (Monte Carlo Markov Chain) algorithm for large V. The number of required 
operations typically grows linearly in V as only sparse matrix type operations are 
required. 

The superfractal idea can be understood as follows. The process of passing from 
the V^-tuple of sets at one level of construction to the y-tuple at the next (Figures 3, 6) 
is given by a random function J-'"' : H^, where is the set of all \^-tuples of 

compact subsets of the plane or of some other compact metric space as appropriate, 
and where a belongs to some index set A. All information necessary to describe the 
chosen J^" at each level in Figure 3 is given by the chosen IFSs for each buffer at that level 
(namely G, F, G, F, F across level 2) and by the arrows pointing to each buffer at that 
level (indicating which three buffers at the previous level are used for each application 
of F or of G), see also the captions in Figure 6. Each has a certain probabihty V"' of 
being chosen, this probability is induced in the natural manner from the probabilities 

and P'-^ of selecting F or G. The JF" are contraction maps on in the Hausdorff 
metric, with contraction ratio equal |, and in general the contraction ratio of J-""- equals 
the maximum of the contraction ratios of the individual maps in the IFSs being used. In 
particular, the superlFS (JF", V , a E A) is an IFS operating not on points in as for 
a standard IFS but on ^-tuples of sets in H^. From IFS theory applied in this setting, 
there is a unique superfractal set and superfractal measure which with probability one 
is effectively given by the collection of ^-tuples of variable fractals together with the 
experimentally obtained probability distribution arising from the previous construction. 
See Barnsley, Hutchinson and Stenfio (2003a) for detailed proofs. 



5. Examples of 2-vetriable fractals 

We begin with 2 IFSs C/ = (/i, /2) ("Up with a refiection") and D ^ {gi, ^2) ("Down"), 
where 

/X 3y 1 X 3y 9 \ fx3y 9 x 3y 17\ 

/i(^,l/)=(2+y-Y^, 2-y+16> ^^^^'^^=l2-T+16'-2-Y+16J' 

fx 3y 1 X 3y 7 \ , , fx 3y 9 x 3y 1 \ 

^^^^'^^=l2+y-16'-2+y+16> ^^^^'^)=l2-T+16' 2+T-16> 

The corresponding fractal attractors are shown in Figure 5. 
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Figure 5. Up (green) and Down (red) attractors 



Figure 6 shows the first 20 steps in the construction of a sequence of pairs of 2- 
variable fractals from the two IFSs U and D. The initial pair of input figures can be 
arbitrarily chosen, here they are each the same and consist of four leaves. 

For the first step in the construction (producing the contents of the second pair of 
buffers) the IFS U = (/i,/2) was chosen, /i was applied to the previous left buffer L 
and /2 was applied to the previous right buffer R; the second buffer was obtained by 
applying U with fi and /2 both acting on the right buffer R at the previous step. Thus 
the first step in the construction can be described by U{L, R) and U(R, R) respectively; 
see the caption below the second pair of screens in Figure 6. The second step is given 
by D{R, R) and D{R, L), the third by U{L, R) and D{R, L), the fourth by U{R, R) and 
D{L, L), and so on from left to right and then down the page. 





Initial Sets 





U{L,R) 



U{R,R) 



D{R, R) 



D{R,L) 
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D{L,R) D{R,R) U{R,L) U{L,R) 




U{L,R) U{L,L) U{R,R) U{R,L) 

Figure 6. A sequence of pairs of (approximate) 2-variable fractals (from left to 
right and then down the page). 

In each case in this particular example, for each buffer at each step, either U or 
D was chosed with probability |. For each buffer, the (in this case two) input buffers 
chosen from the previously generated pair of buffers were also each chosen as either L or 
R with probabihty |, and the same buffer is allowed to be selected twice. After about 
12 iterations, the images obtained are independent of the initial images up to screen 
resolution. After this stage the images (or "necklaces") can be considered as examples 
of 2- variable fractals corresponding to the family [U, D) of IFSs with associated choice 
probabilities (|, |). 

The pair of 2-variable fractals obtained at each step depends on the previous choices 
of IPS and input buffers, and will vary from one experimental run to another. However, 
over any sufficiently long experimental run, the empirically obtained distribution on 
pairs of 2-variable fractals will (up to any prescribed resolution) be the same with 
probability one. This follows from ergodic theory and the fact that the construction 
process corresponds to the chaos game for an IPS (operating here on pairs of images 
rather than on single points as does a standard IPS). As discussed in the previous 
section, we call this type of IPS a superlFS. The collection of 2-variable necklaces 
obtained over a long experimental run should be thought of as a single superfractal, 
and the corresponding probability distribution on necklaces should be thought of as the 
corresponding superfractal measure. 

In Pigure 7 we have superimposed the members of a generated sequence of 2- 
variable fractal necklaces. By virtue of the fact that, as discussed before, the probability 
distribution given by such a sequence approximates the associated superfractal measure, 
the image can be regarded as a projection of the superfractal onto 2-dimensional space. 
The attractors of the individual IPSs U and D are shown in green and red respectively. 
The projected support of the superfractal is shown on a black background but, inside 
the support, increasing density of the superfractal measure is indicated by increasing 
intensity of white. 
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Figure 7. Superfractal projected onto 2 dimensional space. 

In Figure 8 are shown a fern and lettuce generated by two IFSs, each IPS consisting 
of 4 functions. In Figure 9 is a sequence of hybrid offspring, extending the examples 
in Figure 4. The colouring was obtained by working with two IFSs in 5 dimensional 
space, with the three additional dimensions corresponding to RGB colouring. The two 
IFSs used project onto two IFSs operating in two dimensional space and which give 
the (standard black and white) fern and lettuce attractors respectively. The 2-variable 
offspring were coloured by extending the superfractal construction to 5 dimensional 
space in a natural manner. 




Figure 8. Lettuce and fern attractors. 
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Figure 9. A sequence of fern-lettuce hybrid offspring. 



6. Computation of Dimension for \^-variable Fractals 



An important theoretical and empirical classification of fractals is via their dimension. 
We first show how to compute the dimensions of \^-variable Sierpinski triangles. As we 
discuss in the next section, the method generalises. 

Associated with the transition from the fc-level set of V buffers to the k+1 level is 
a matrix M^{a) defined for each a as follows. Entries of M''{a) are initialised to zero. 
For the set in the f th output buffer at level k + 1 one considers each input buffer w used 
in its construction and adds to the wth entry in the vth row where r = | or | is 
the corresponding contraction ratio. The construction of M''{a) can be seen in passing 
from level 1 to level 2 and from level 2 to level 3 in Figure 3, giving respectively: 
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The "pressure" function 

7v(a) = ^lim ^ log (^^ \\M\a) • . . . • M'{a)\\^ (2) 

exists and is independent of the experimental run with probability one by a result of 
Furstenberg and Kesten (1960), see also Cohen (1988) for the version required here. 
(By ||74|| we mean the sum of the absolute values of all entries in the matrix A.) The 
factor 1/y is not necessary in the limit, but is the correct theoretical and numerical 
normalisation, as we see in the next section. (See Feng and Lau (2002) for another use 
of Furstenberg and Kesten type results for computing dimensions of random fractals.) 
In case V — 1, 

7i(«)=(l-f)log3-|log2 (3) 

from the strong law of large numbers. 

It can be shown, see Barnsley, Hutchinson and Stenflo (2003b), that for each V 
7v(q;) is monotone decreasing. In this example the derivative lies between — log 2 and 
— log 3, corresponding to the contraction ratios | and | respectively, see Figure 10. 
Moreover, there is a unique d = d{V) such that 7y((i) = 0. This is the dimension of the 
corresponding V"-variable random fractals with probability one. The establishment and 
generalisation of this method uses the theory of products of random matrices and ideas 
from statistical mechanics, as we discuss in the next section. 
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Figure 10. Graphs of the "pressure" function 7y(a) for V = 1, 2, and 5 
respectively, from left to right. 
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It was previously known that the dimension of homogeneous random Sierpinski 
triangles {V = 1) is 21og3/(log2 + log 3) ^ 1.226, see Hambly (1992, 2000), and that 
the dimension of standard random Sierpinski triangles {V — > oo) is the solution d of 
^3 (I) + |3 (I) = 1, or approximately 1.262, see Falconer (1986), Graf (1987) and 
Mauldin and Williams (1986). In particular, (jH} is in agreement when computing d{l). 
For \^ > 1 we used Monte Carlo simulations to compute 7y(a) in the region of the 
interval [d{l),d{oo)], with the computed values shown (Figure 10). These values have 
error at most .001 at the 95% confidence level, and from this one obtains the dimensions 
d{2) ^ 1.241, d{5) ~ 1.252 (Figure 10). The computed graphs for V > 1 are concave 
up, although this does not show on the scale of Figure 10. 

7. Analysis of Dimension results for ^-variable Fractals 

In order to motivate the following analysis, consider a smooth curve or smooth surface 
having dimension 1 or 2 respectively. It is possible to cover each "efficiently" (i.e. with 
little overlap) by sets of small diameter such that the sum of the diameters raised to 
the power 1 or 2 respectively is very close to the length, or to the area divided by 7r/4, 
respectively. However, if any power a > 1 or 2 respectively is used then the limit of this 
sum, as the maximum diameter of the covering sets becomes arbitrarily small, is zero. 
For any power a < 1 or 2 respectively the limit of the sum, as the maximum diameter 
becomes arbitrarily small, is infinity. 

In the case of the construction of 5-variable Sierpinski fractal triangles, we see that 
if one begins the construction process with 5 copies of a triangle T as indicated (Figure 
3), then the contents of each buffer at later stages will consist of a large number of 
tiny triangles which approximate and cover the "ideal" or limiting 5-variable Sierpinski 
fractal triangles. (Note that what is obtained after k steps is only an actual 5-variable 
Sierpinski fractal triangle up to k levels of magnification.) One can check that the sum 
S^{a) of the diameters to the power a of the triangles in the v-th buffer at level k is 
given by the sum of the entries in the v-th row of the matrix M^{a) ■ . . . ■ M'^{a). 

Using (j21) it is not too difficult to show that lim^^oo l ^og S^{a) also exists for each 
V and equals 7y(a), independently of v, with probability one. (The argument relies on 
the existence of "necks" , where a neck in the construction process occurs at some level 
if the same IFS and the same single fixed input buffer is used for constructing the set 
in each buffer at that level.) One can show (Barnsley, Hutchinson and Stenflo 2003b) 
that 7y(Q;) is decreasing in a and deduce that there is a unique d = d(y) such that 
Jvid) = (Figure 10). From this and Q it follows that 



It is hence plausible from the previous discussion of curves and surfaces that the 
dimension of V^-variable Sierpinski triangles equals d{V) with probability one. The 
motivation is that the covering by small triangles is very "efficient" . 
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The justification that the dimension is at most d is in fact now straightforward 
from the definition of (Hausdorff) dimension of a set. The rigorous argument that 
the dimension is at least d, and hence exactly d, is much more difficult, see Barnsley, 
Hutchinson and Stenflo (2003b). It requires a careful analysis of the frequency of 
occurrence of necks and the construction of Gibbs type measures on variable Sierpinski 
triangles, analogous to ideas in statistical mechanics. 

Similar results on dimension have been established much more generally, see 
Barnsley, Hutchinson and Stenflo (2003b). For example, suppose the functions in each 
IFS are similitudes, i.e. built up from translations, rotations, reflections in lines, and a 
single contraction around a fixed point by a fixed ratio r (both the point and the ratio r 
may depend on the function in question). We also require that the IFSs involved satisfy 
the uniform open set condition. In the case of 5-variable Sierpinski fractal triangles 
constructed from the IFSs F = (/i, /2, /s) and G = {gi, g2, gs) this means the following. 
There is an open set O (the interior of the triangle T) such that fi{0) C O, f2{0) C O, 
fz{0) C O, and fi{0), f2{0), f^iO) have no points in common, and such that analogous 
conditions apply to the maps gi, g2, gs with the same set O. Under these circumstances 
one constructs the matrices M^{a) and the pressure function 7^(0;) as before and it 
follows that the solution d{V) of 7y(d) = is the dimension of the corresponding V- 
variable fractals with probability one. 

8. Generalisations 

Many generahsations are possible. 

The number of functions M in each of the IFSs F" = (/", . . . , . . . ,Pm) 

vary. The maps need only be mean contractive when their contraction ratios are 
averaged over m and n according to the probabilities and P"^ respectively. Neither the 
number of functions in an IFS nor the number of IFSs need be finite; this is important 
for simulating various selfsimilar processes, including Brownian motion, sec Hutchinson 
and Riischcndorf (2000). The maps may be nonlinear, and many of the results and 
arguments, including those concerning dimension, will still be valid. The set maps 
need not be induced from point maps; this is technically useful in extending results to 
the case where M is not constant by artificially adding set maps such that f^{A) is 
always the empty set. It could also be important in applications to modelling biological 
or physical phenomena where the objects under consideration are not naturally modelled 
as sets or measures. Buffer sampling need not be uniform; buffers could be placed in a 
rectangular or other grid, and nearby buffers sampled with greater probability, in order 
to simulate various biological and physical phenomena. 

An IFS operates on R^, or more generally on a compact metric space (X, rf), to 
produce a fractal set attractor; a weighted IFS produces a fractal measure attractor. 
We have seen in this paper how a family of IFSs operating on (X, d), a probability 
distribution on this family of IFSs, and an integer parameter V, can be used to generate 
a (super)IFS operating in a natural way on (IHI(X)^, ^7^), where IHI(X) is the space 
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of compact subsets of X and dy^ is the induced Hausdorff metric. In the case of a 
family of weighted IFSs, the induced supcrlFS operates on (F{X)^ ,dMK) where P(X) is 
the space of unit mass measures on X and dMK is the induced Monge Kantorovitch 
metric. In either case there is a superfractal set (consisting of V^-variablc sets or 
measures respectively) together with an associated superfractal measure (a probability 
distribution on the collection of variable sets or measures) . There is also a fast forward 
algorithm to generate this superfractal. 

We can consider iterating this procedure. Replace (X, d) by (]HI(X)^, dn) or 
(P(K)^ ,dMK): take a family of superlFSs, an associated probability distribution, and 
a new parameter W. To speculate; if a superfractal may be thought of as a gallery of 
a new class of fractal images, can one use some version of the iteration to produce a 
museum of galleries of yet another new class of fractal images? 

9. Conclusion 

There appear to be many potential applications, which include both the extension of 
modelling possibilities to allow a controlled degree of variability where deterministic or 
random fractals have been previously applied, and the rapid generation of accurately 
distributed examples of random fractals — previously not possible except in very special 
cases. 
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